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Abstract 

We consider wave propagation and scattering governed by one dimensional Schrodinger 
operators with truncated periodic potentials. The propagation of wave packets with 
narrow frequency supports is studied. The goal is to describe potentials for which the 
group velocity (for periodic problem) is small and the transmission coefficient for the 
truncated potential is not too small, i.e. to find media where a slowing down of the 
wave packets coexists with a transparency. 
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I. Introduction 

This paper is designed to provide a mathematical foundation for studying the processes 
of slowing down of the wave packets in periodic media. The problem of creating compact 
and efficient optical delay devices has been very actively discussed in literature for the last 
several years, but it has not yet been solved technologically. Such devices can find many 
applications, for example in synchronizing the work of very fast optical elements and much 
slower electronics. Several different mechanisms to slow down the propagation of wave pack¬ 
ets have been considered (see iia, lu, n, i, a, a, a, ai). Most of them are based 
on a use of periodic media with a flat dispersion relation lo = uj{k) in some region. For such 
media, the derivative k'{u!) is very big, and the group velocity Vg{u!) = is small for a 
specihc interval of frequencies \u} — uJo\ < e. Here L is the period. 

Let iL be a periodic 1-D Schrodinger operator on L‘^{R) : 

(R'ip 

H'tp = — r^ + q{x)ip, q{x + L)=q{x). (1) 

dx^ 

Let b = be a Bloch solution of the equation Hip = oj'^ip, i.e. ip{x + L,uj) = 

ip{x^ijj)e''^^^\ The function k = k{uj) defines the dispersion relation, and are the 

eigenvalues of the monodromy operator ilU. The function ip{x,uj)e~^^^ is a solution of the 

* Corresponding author. E-mail address: brvainbe@uncc.edu (B. Vainberg). Both authors were partially 
supported by the NSF grant DMS-0405927. 


1 





nonstationary wave equation uu = —Hu, it corresponds to a wave which propagates along 
the x-axis with the (phase) velocity V{uj) = uiL/kiui). Consider a wave packet 


u = 


ip{x,uj)a{uj — uiQ)e 


\u)-Wo\<£ 


( 2 ) 


with angular frequencies from a small interval around (Jq. The group velocity for the packet 
m is defined a.s Vg = L/k'^uJo). One can easily show that 


where a is the Fourier transform of the function a (a;), 7 is a periodic function, i.e. u has 
the form of a plane wave which propagates with the phase velocity V{ijJq) and which is 
modulated by another wave 5 whose speed of propagation is Vg. Asymptotic formula © is 
valid when 5 and |a:| are small. The main physical idea behind the optical delay devices 
is to find an optical material for which the group velocity Vg is very small on some interval 
around a specific ’’singular” frequency ojq and then to use the wave 5 (rather than the last 
factor in the right hand side of ©) for the transmission of information. 

It is well known that the spectrum of the periodic operator H has a band-gap struc¬ 
ture with the bands defined by the equation |Trace M^\ < 2. The function k{uj) changes 
monotonically on any band, the changes are equal to ±7r, and therefore, k' is large and Vg 
is small on any narrow band. This suggests solving the problem of slowing down of the 
wave packets by creating a periodic optical material with narrow spectral bands. A group 
of physicists from Caltech (A. Yariv and others [1^1, [H], [I2I) proposed to use a special 
coupled resonators optical waveguide (CROW) as an optical delay device. It can be realized 
as a set of identical periodic cavities (arranged in a linear sequence) in a dielectric medium. 
If there are no interactions between resonators, then the spectrum of the system consists of 
discrete eigenvalues (of the problem in an individual cavity) of infinite multiplicity. If a weak 
and periodic interaction between neighboring resonators is imposed, then such eigenvalues 
will generate narrow optical bands. P. Kuchment and A. Kunyansky [U] provided spectral 
analysis of very general high contrast periodic media which yielded many examples of high 
contrast media with narrow spectral bands (and therefore, with the flat dispersion relations 
and small group velocity). 

Hence, one can expect that periodic media with small group velocity (for example, due to 
the existence of a narrow spectral band) may serve as optical delay devices. However, let us 
note that all the arguments above are related to infinite periodic media, and a real device has 
finite size. The main goal of this paper is to answer the following question: in what sense may 
the problem in an infinite periodic medium serve as an approximation for the processes of 
wave propagation through a long finite slab of the periodic medium. This question becomes 
particularly important due to the following circumstance: one needs to know if the incident 
pulses with the frequencies supported in a small neighborhood of a singular frequency ujq 
(where the dispersion is flat) will enter the slab or, perhaps, the majority of their energy will 
be reflected. 

In this paper we study a Hamiltonian of the form 


Hn^^ 


dx‘^ 


qN{.x)tlj, 


(4) 
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where the potential coincides with a periodic real valued potential q{x) when x E [0, iVL], 
and qjv = 0 outside the interval [0,iVL], and N is very big. The goal is to hnd a potential 
q such that for some hxed interval A on the frequency axis the group velocity for operator 
m is small and the transmission coefficient for operator m is not too small. We shall see 
that in order to achieve this goal we need to consider a family of operators with potentials 
q{x) = Av{x), A —>■ cxo. The parameter A has to be large in order to make the group velocity 
small. At the same time the interval A shrinks as A —>■ cxo, and therefore, we need to assume 
that A —> cx) (a physical motivation for that will be given). We shall see that there are two 
scenarios with similar results: A is a narrow band for the operator m or it is a neighborhood 
of a degenerated band edge. The dispersion relation is very flat in the hrst case, and it is 
not flat in the second case. In both cases, the length of A and the group velocity have order 
0{A~^). The transmission coefficient has the same order on the main part of the interval A, 
but there are N very narrow transparency zones on A. We shall also justify the fact that in 
spite of these transparency zones, only a small part of the energy (of order 0{A~^)) of the 
incident pulse will be transmitted through the media, and the majority of the energy will be 
reflected. 

We shall not discuss other known models of slowdown of wave packets: inflection point 
(A. Figotin, I. Vitebskii, j2]), necklace model (S. Molchanov, B. Vainberg P), SCISSOR (J. 
Heebner, R. Boyd, i.ffl) mainly because the corresponding mathematical tools are different 
from the classical 1-D scalar Schrodinger operator: it is symplectic systems in P, quantum 
graphs in P, etc. 

The paper has the following structure. 

In the second section, we develop the scattering theory for the operator (ED combining 
the ideas of the Floquet-Bloch theory and the 1-D scattering theory for operators with 
fast decaying potentials. In particular, exact formulas for the reflection and transmission 
coefficients will be obtained through the monodromy operator and N. Section II contains 
also a transition to semi-inhnite periodic materials {N —> cxo). The relation between a 
slowing down and the reflection is discussed in the last section. 

II. Scattering by a finite slab of a periodic medium. 

a). General Schrodinger operators and operators with periodic or compactly supported 
potentials. We shall hrst recall some well known facts concerning the spectral problem for 
the general Schrodinger operators on {R) : 

Hxjj = —ijj'' -\- q{x)ip = (5) 

where the potential q (which can be a distribution) is real valued and bounded from below 
in an appropriate sense, say, there is a constant c such that 

x+l 

min(g (x), Qi)dx > c for all x E R. (6) 

The last inequality provides the uniqueness of the self adjoint extension in (i?) of the op¬ 
erator H dehned originally on the space of inhnitely smooth, compactly supported functions 
ijj. 
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Let (0, x) be the Priiffer transfer matrix (propagator) for operator H 


T (a ^ 

where 2 i^y the solutions of the equation Htl^ = uj'^'ip with initial data 

V’i(o^,0) = l, ij[{uj,0) = 0; V'2(^,0) = 0, i/j^iuJ^O) = I, 


(7) 

( 8 ) 


(i.e., simply, (0, 0) is the identity matrix). For any solution -0 of the equation = a;^'0, 
matrix (0, x) maps the Priiffer Cauchy data of -0 at a; = 0 into the Priiffer Cauchy data 
of ip at point x : 

p) [x) 

LO 

Assumption © implies that the equation Ilpj = oj‘^ip with Imcn > 0 has exactly one 
solution in {R+) normalized by the condition [u, 0) = 1, and it has exactly one 
solution 'tp~ G {R-) normalized by the same condition. Here R± are the semiaxes x ^ 
0. Obviously, can be represented as linear combinations of V’l and ip 2 y and from the 

normalization of ^jJ^ it follows that there exist functions (a;) such that 


7b (0, x) 


b(0) 

h'(o) 


ijj^ = bi + (<^) b’2) > 0. 


In other words, are functions such that 


bi + rnR (a;) -02 ^ {R+) y + m (u) ^2 ^ {R-) y Icncn > 0. 

Functions (a;) are called Weyl’s functions, and are called Weyl’s solutions. 

Let the potential q = q{x) be periodic: q{x + L) = q{x). Consider the propagator through 
one period (monodromy matrix): 


(0, L) 


a (3 

7 5 


(^) 


Put F {ui) = -TrM^ = + ip' 2 ) {ix^L) = (a + 5)(ci;). Both and F (uj) are entire 

functions of oj. Since det = 1, the eigenvalues /i^ (cn) of are the roots of the following 
characteristic equation: 

/i^ + 2F (w) /i + 1 = 0. (9) 

If Irncn > 0 then one can select roots /i^ (a;) of m in such a way that /i^ (a;) = 
where k (a;) is analytic and 


i.e.. 


Imfc (cj) > 0 when Imcn > 0, 


l/i’*' (a;)I < 1, |/i (o')I > 1, Imcn > 0. 


( 10 ) 

( 11 ) 


The roots (u) for real a; > 0 are dehned by continuity in the upper half plane: 

/i^ (u) = (u + iO), u E [0, cx)]. 
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Figure 1: A typical graph of F (uj) . Segments bi = [a;i,a; 2 ], &2 = [<^ 3 , <^ 3 ], ^3 = [ 1 ^ 5 , i^e], 
64 = [cuy, Co's] are bands, the point ojq = (Uy is a degenerate band edge. 

Since the trace of is equal to the sum of the eigenvalues 

cosfc(a;) = F (a;) = (a + (5)/2. ( 12 ) 

If the potential q is not negative then the spectrum of H belongs to the positive part of 
the energy axis A = o;^, and the spectrum has a band-gap structure. If g_ = min(g, 0) 7 ^ 0 
then H may also have spectrum (a finite number of gaps, in particular) on the negative part 
of the energy axis A = —o;^ < 0. Wave processes are associated only with the positive part 
of the spectrum, and we shall exclude the negative part of the spectrum from the future 
analysis. Thus, the bands on the frequency axis uj will always be related to the positive 
energies. 

For real a; > 0, the inequality |F(a;)| < 1 defines the spectral bands (zones) = 
[a; 2 n-i, <^ 2 n ], n = 1,2,..., on the frequency axis a; = y/X. The function k{uj) is real valued 
when UJ belongs to a band. The roots (a;) are complex adjoint there, and |/i^| = 1. The 

00 

spectrum of H (on (i?)) on the frequency axis is IJ bn- The complimentary open set, 

n=l 

given by \F (a;)| > 1, corresponds to spectral gaps. On gaps, the roots /i^ (a;) are real and 
m holds. Figure 1 presents a typical graph of F (a;). A point ujj which belongs to the 
boundary of a band and the boundary of a gap is called a non-degenerate band edge. If it 
belongs to the boundary of two different bands, it is called a degenerate band edge. 

Lemma 1 If uj = ujq is a non-degenerate band edge, then F'{uJo) 0. If uj = ujq is a 

degenerate band edge, then F'{uJo) = 0, F"(a;o) 7 ^ 0. Both eigenvalues of the monodromy 
matrix at any band edge are equal to 1 or both are equal to —1. The matrix at the 
non-degenerate band edges has an off-diagonal element in its Jordan form, and = ±J at 
the degenerate band edges, where I is the identity matrix. 

We shall provide only a sketch of a proof of this Lemma, and we refer to [Hj for details. 
Let us recall that bands are defined by the condition |F(a;) | < 1, and that F (uj) = ±1 
at any band edge u = uq- If F'{uJo) 7 ^ 0 then |F(ci;) | — 1 changes sign at uj = ujq, and 
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a; = a;o is a non-degenerate band edge. Let a; = a;o be a band edge and F'{ujq) = 0. Using 
the equation in variations, one can show that in this case F”{u!q) ^ 0 and F{ujo)F"{u!o) < 0 . 
Hence, \F (a;) | < 1 in a neighborhood of such a band edge, i.e. a; = cuo is a degenerate band 
edge. 

The eigenvalues of are equal to where cos k (a;) = F (a;) = ±1 at a band 

edge uj = ujq. Hence, both eigenvalues at a; = ojq are equal to 1 or both are equal to 
— 1 (and F{(jJq) = 1 or F{ijJq) = —1, respectively). Thus, A = cUq is an eigenvalue of the 
periodic or, respectively, anti-periodic problem for the operator H. If the band edge is non¬ 
degenerate, then this eigenvalue is simple, and the Jordan form of TU, is non-diagonal (since 
the eigenspace is one dimensional). If the band edge is degenerate, then the corresponding 
eigenspace is two dimensional, and the Jordan form of is diagonal, and therefore, TU, = 
±. 


The proof of Lemma ^ is complete. 

We normalize the eigenvectors (a;) of by choosing the first coordinate of (ci;) 


to be equal to one: (cj) = 


(u) 


The second coordinates of the vectors (oj) 


coincide with Weyl’s functions dehned above. In fact, if "0^ are solutions of the equation 
with the Priiffer Cauchy data given by the eigenvector h^, then 


0^ (a;, x + L)= x ), 


(13) 


and CO implies that 0^ G {R±) when Imuj > 0. From here it follows that 0^ coincide with 
Weyl’s solutions introduced for general Hamiltonians iJ, and that the second coordinates of 
the vectors (cu) are Weyl’s functions. 

Since 


■ a - 0 


■ 1 



mA 


(14) 


the following two representations are valid for Weyl’s functions: 


m 


[UJ] = 


g±ifc(a;) _ ^ 


7M 


- J (a;)' 


(15) 


For real a; G IJ f^e Weil’s solutions 0^ (a;,a;) are the usual Bloch almost periodic 

n=l 

oo 

solutions. The relationship between the space and time frequencies k = k (a;), a; G u bn, is 

n=l 

called the dispersion law. 

Consider now a Hamiltonian H with a compactly supported potential q. Let suppg C 
[0,H]. A function 0 is called the solution of the scattering problem for the operator H if 
iL0 = and 

0 = 0"’^ + re-®^0 x<0; = x > A, 

with some r = r(a;) and t = t{u}). The numbers r and t are called the reflection and 
transmission coefficients, respectively. 

For any matrix M = let ||M|| be the Hilbert-Schmidt norm of M : 


||M|p = 
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Theorem 2 Let the support of q belong to [0, A]. Then the following relations hold for the 
reflection and transmission coefficients: 




\\TU0, A)\\^ - 2 

|| T .( 0 , A )||2 + 2 ’ 




4 

\\tuo,aw + 2- 


(16) 


Proof. The Priiffer Cauchy data for the scattering solution at x = 0 and x = A are 


1 + r 


t 

_ z (1 - r) _ 


it 


respectively. This implies 


t 

it 


(0, A) 


1 + r 
z (1 - r) 


We equate the ratios of the second and the hrst coordinates of the vectors above. Using 
the notation 

a b 


U.(0,/1) = 


c d 


(cn) 


we arrive at 


c(l + r) + di{l — r) 
a(l + r) + bi{l — r) 


Solving for r we obtain 


r = 


Thus, 


\rf = 


{d 


d — a — i{c + d) 
a + d + i{c — d) 

(17) 

+ dy \\T^ {0, A) y — 2ad + 2cd 
— dy \\Ti^ {0, A) y + 2ad — 2cd 

(18) 


Since detT^ (0,^) = 1, implies the first of relations (ITHll . The second relation follows 
from the first and the energy conservation law: |rp + \t\^ = 1. 

The proof of the Theorem is complete. 

b) Scattering by truncated periodic potentials. The main result of this section concerns 
the Hamiltonian with the truncated periodic potential qn '■ 


HN'f = 


d^ 

dx'^ 


+ qN{x)'f, qN{x) = q{x)fo^NL], 


which appears when the propagation of waves through a hnite slab of a periodic medium is 
studied. We shall also consider the limiting case N = oo : 


Hoof 


dx'^ 


+ qoo{.x)tlj, 


qoo{x) = g(x)/[0,oo], 


which corresponds to the case of such a long slab that it can be considered as half infinite. 

We shall denote the reflection and transmission coefficients for the operator Hn (with a 
compactly supported potential qisf) by rjv and t^, respectively. In the case of the operator 
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Hoo, the solution ifj of the scattering problem is defined as the solution of the equation 
Hoo'4^ = which has the form 

if) = a; < 0; '?/> = x), x > 0, (19) 

with some r = r(a;), c = c(lij). 


Theorem 3 Ij. The transfer matrix over N periods Mff = (0, LN) has the form 

sin Nk (cn) sin {N — l)k (cn) 

sin k (cn) sin k (cn) 

where k = k (cn) is the dispersion relation. The elements of Mff satisfy the relations 


= 


fdjyl 

In 


( 20 ) 


C^N — Sn 


In 


sin Nk (cn) 
sin k {oo) 
sin Nk (cn) 
sin k (cn) 


(tt - 5 ), I 3 ^ = 

7 , otv + = 2 


sin Nk (cn) 
sin k (cn) 


/9, 


cos Nk{(jj). 


( 21 ) 


The Hilbert-Schmidt norm of Mff is equal to 


m:: 


(iiMjr 


sin^iVfc^ 
^ sin^k{u) 


2). The reflection coefficients have the forms 


tn (c-^) 


r (cn) 


{a - S) + i {P + ■j) 

2sinfc(o;) ZTkt) + i - PY 

/3 + 7 — i (a — (5) 

2 sink (cn) + (/5 — 7 ) ’ 


3). The following formula is valid for \tN {eo) 


\tN (oj) 




4 

r)\sin^ Nk(u>) . a 
sin^ fc(a;) + ^ 


( 22 ) 


(23) 

(24) 


(25) 


4). The reflection coefficients r^ (cn) and r (cn) are analytic in the upper half plane C~^ = 
{u :Imk > 0} and continuous in C+, and un (oj) —^ r (u) when N ^ 00 and u E C^. When 
u is real, r^ (cn) converges to r (cn) in the weak sense: 


rjsf {oj) ip{(jj)duj 


r (cn) (p(u)duj as N ^ 00 , 


for any test function ip E D. 



Remark. The convergence of rN{uj) to r (a;) follows from explicit formulas but 
it also can be derived from a much more general fact that the spectral measure of the 
Schrodinger operator H in L?‘{R) is a weak limit of the spectral measures of the operators with 
truncated potentials. This fact is valid for very general potentials without the assumption 
on periodicity. We shall discuss this result elsewhere. 

Proof. The monodromy operator satishes the relation Q: 

— 2 cos k (a;) + / = 0 . 

Formula dm) follows from here by induction. The hrst three relations of m are immediate 
consequences of P|). In order to get the fourth one we note that the eigenvalues of are 
(u) = Thus, 

aN + SN = TrM^ = = 2 cos Nk (uj). (26) 

In order to prove (1^ we note that the relations det = det = 1 imply that 

||MJp-2 = (a-^)2 + (/3 + 7)^ (27) 

and a similar relation is valid for ||iW^|p — 2. From here and (I21|l it follows that 

lKir-2 


This justihes (1^ . 

Let us prove the second statement of the Theorem. From (ED and m it follows that 

+ + ((a - ^) + »()3 + t)) 

^ + Sn + t (7n -2cosNk h - P) 

This justihes (El- In order to get (1^ we note that (USD implies that 


(«^ - + (/3^ + -S) + {P + 7)^] 


sin^ k (u) 


sin^ Nk (o;) 
sin^ k {oo) 


[iiMjr-2]. 


or 


i(l-r) 


1 + r 


1 

_ i [1 - r] _ 

= c 

m'^ 


and therefore, r = \jJ^+ ■ From here and |T5|l it follows that 


r 


r 


i — 
i + 
i — 
i + 


0 

/? 


7 

7 

e^^ — 5 


a + i(3 




— {a — i(3) ’ 

— 6i — j 

i^lk _ _j_ ry 


and 

— (5 — Z7) 
g*A: — (5 + ^y) 


Hence, 

^ ^ (q + i(5) -{5- i-f) ^ (a - 7 + i (/3 + 7 ) ^ (a - ^) + i (/3 + 7 ) 

2 e*^ — (a — iP) — {6 + ij) 2 e*^ — (a + 5) + i (/3 — 7 ) 2ismk + i {P — 7 ) 
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The last equality is a consequence of m, and it implies (12411 . 

The third statement of the Theorem follows immediately from (ITT)|l and (j22I)- The ana- 
lyticity of tn {^) and r in and their continuity in C+ follow from the explicit formulas 
(ESI), (El- Furthermore, if Imk (a;) > 0, then 


cos Nkiuj) 
sin Nk{ui) 


as N 


CX), 


and this justihes the convergence oir^ (cc>) to r when u E The weak convergence on the 
real axis is an obvious consequence of the convergence in the complex half plane. 

The proof of Theorem El is complete. 

III. Relation between slowing down and reflection. We start this section with 
considering the scattering problem for the operator II]y. The group velocity (for the periodic 
in R potential q) is dehned (see (E))) as 


y, = L/k' (u ). 


(29) 


In spite of the fact that formula Ei concerns periodic in R potentials, one can expect that 
a hnite slab of a periodic medium with a small group velocity, defined by (ESI), can be used 
as a device for slowing down the propagation of wave packets. However, it could happen 
that the wave packets with frequencies around u = Ug, for which Vg (ug) is very small, will 
not enter the slab, but will be reflected almost completely. Hence, one needs to And the 
situations when simultaneously Vg is small and the reflection is not very large. To study this 
problem we consider the solutions of the stationary scattering problem for the operator Rjv 
with values of a; = ug, for which Vg is small, and evaluate the reflection coefficient rjv(u’g). 
When it is needed, we shall make more rigorous analysis of the situation by studying the 
corresponding time dependent problem. 

Let us discuss the reasons for the group velocity to be small. Due to m, the dispersion 
k = k (o;) on each band bn is equal to ± arccosF (a;), where F (a;) = ^TrM^. The sign in 
this representation depends on the sign of F' (oj). Hence, 


k' (u) 


±F' (u) 


(30) 


and k' (a;) can be large either because oj is close to a value oog for which |F(a;o) | = 1 or 
because |F’'(a;)| >> 1. In the first case, there are still two different possibilities: F' {oog) ^ 0 
or F' (ojg) = 0. If F' (ojg) 7 ^ 0 then oog is a non-degenerate band edge (see Lemma [T)). If 
F' (ojg) = 0 then oog is a degenerate band edge: two bands are adjacent at oog, but an 
arbitrarily small perturbation of the operator may open a gap around ojg. One can consider ojg 
as a degenerate gap consisting of one point (see Fig. 1). Thus, we should study the following 
three cases: the spectrum of the incident wave packet is supported in a neighborhood of a 
non-degenerate band edge, in a neighborhood of a degenerate band edge or in a region where 
|F'(a;)| » 1. 

Wave packets with frequencies in a gap or near a non-degenerate band edge. Note that 
the transmission coefficient tv(<^) can not vanish at a fixed a; if < oo, since otherwise the 
scattering solution is zero for x > NL, and therefore has to be identically zero. A part of 
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the energy of an arbitrary incident wave packet will be transmitted through a hnite slab of 
any periodic medium. However, the following statement shows that the transmission will be 
negligibly small if the slab is long [N is big) and the frequency spectrum of the incident wave 
packet is supported in a gap of the corresponding periodic problem or near a non-degenerate 
band edge. Let us recall that if ujq is a band edge, then | cos/c(a;o)| = 1, i.e. k{uJo) = nn, n 
is an integer. 


Theorem 4 1). If ut belongs to a gap of the periodic Hamiltonian (0)) then 

|fAr(a;)| = as N ^ oo, 


where a > 0 depends on the distance from uj to the closest band. 

2). IfujQ is a non-degenerate band edge andoj is so close to Uq that N\k{uj) — k{u!o)\ < 7r/2 
then 


l^vMI = 0(—) as N 


oo. 


(31) 


Remark. The assumption in the second statement means, roughly speaking, that \uj — 
a;o| = 0(l/iV). 

Proof. If UJ belongs to a gap then pH = are real, i.e. k (a;) is purely imaginary 

modulus nn. Thus, 


I sin Nk (a;) | 

I sin k (a;) | 

If ujq is a band edge, then |/i=*"(a;o)| = 


= 0(e'"^) asN ^oo. 

1 and k (a;o) = nvr. If N\k (a;) 


(32) 


nvrl < 7r/2 then 


sm. Nk (uj) 
I sin k (a;) | 


0(N) as N ^ oo. 


(33) 


The statements of the Theorem will follow from (IS21), (ISSI), and if we show that 
I |Mt^| p > 2 when uj belongs to a gap or a neighborhood of a non-degenerate band edge. 

Formula (EID implies that for any uj we have > 2, and = 2 only if a = 5, 

f3 = —7. In the latter case, the eigenvalues pH of the matrix are equal to a ± if3. If 
UJ belongs to a gap, then p^ have to be real and (ITT| has to hold. This is impossible if 

= a ± ip. Hence, > 2 in gaps. If ujq is a band edge, then |p=*"(a;o)| = 1. Thus, 

a = 6 = ± 1 , P = —7 = 0 at a; = ci^o. This contradicts Lemma ^ in the case of a non¬ 
degenerate band edge. Thus > 2 and therefore this is also true for uo close enough 

to Ci^o- 

This completes the proof. 

Wave packets with freguencies near a degenerate band edge. We shall show that the 
reflection coefficient r^ is zero at any degenerate band edge uo = uoq, i.e. the medium is 
transparent for the plane wave with the frequency uoo (not for the wave packets). It follows 
immediately from dSOD and Lemma C] that the group velocity is zero at non-degenerate band 
edges. This is not true for degenerate band edges. We shall show that the group velocity 
in this case is ’’usually” not small, but can be made small for a specihc medium. The next 
theorem provides formulas for r^ and Vg at uj = ujq. It is followed by a discussion of how 
the group velocity can be small, and then by a discussion of the transparency of the media 
for wave packets. 
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Theorem 5 Let u = uq be a degenerate band edge (i.e. F (a;o) = ±1, F' (a;o) = 0 Then 


r„(c.„) = 0. and 

Proof. If ci; = a;o is a degenerate band edge then (see Lemma H)) = ±I and 

F" (cuq) 7 ^ 0, where I is the identity matrix. The hrst relation allows us to pass to the limit 
in (El as a; —*• ujq. The numerator in the right hand side of El vanishes as a; —>• ujq. 
The denominator converges to -F2/N, since kiujo) = mi. This proves the hrst statement of 
Theorem El Let us evaluate the group velocity. We shall assume that F (a;) = 1. The other 
case {F (a;) = — 1) can be treated similarly. Then we have: 

F (cn) = 1 — c{ijj — ujof + O (Icn — cuol^) as a; —> cuo, where c > 0. 

This and El imply that k (cn) = 2mi ± (cn — + 0{\ui — where the sign has 

to be chosen to satisfy El- Hence, 

k (a;) = 2mi + (a; — a;o)\/^ + 0(|a; — a;o|^^^), 

and k' {oj) = \f^c + 0 {\uj - = ^/\F%J^\ + 0{\uj — This completes the proof 

of the Theorem. 

Since we are looking for media with a small group velocity, we need the value of F"{ujo) 
to be large. Hence, the potential has to be large, since the elements of the monodromy 
matrix and their derivatives with respect to u can be estimated through the L^-norm of the 
potential over the period. We restrict ourselves to considering a natural class of potentials 
of the form q{x) = Av{x), where H > 0 is a big parameter and v{x) is a piece-wise inhnitely 
smooth periodic function which is separated from zero on each interval of the continuity. We 
shall show that the group velocity is not small if v{x) is positive, negative, or changes sign 
one time on the period. However, a more sophisticated choice of v may provide an example 
of a potential for which |F"(a;o)| —> cxo as H —>■ oo. We shall provide one such an example, 
in which the intervals where v is positive or negative are separated by the intervals where 
n = 0. To make the calculation simpler we consider a non-smooth v of the form 

OO 

q{x) = A [(5(x — n/) — (5(x — 2n/)], / = L/2, (34) 

n=—oo 

where 5(a;) is the Heviside delta function. Let us also mention that the wave equation 
Utt + Hu = 0 is normalized in such a way that the speed of light is equal to one. So, when 
we speak about a small group velocity, we mean that it is small compared to one. 

Theorem 6 Let the frequency oj = ojq be fixed and let the potential q have the form q{x) = 
Av{x), where v is a piece wise C°°function, v{x + L) = v{x), and H —> cx). Then 

a) if V & C°° on (0, L) and v{x) > c > 0, then ojq belongs to a gap when A is large 
enough; 

b) if V E C°° on (0, L), v{x) < —c < 0, and ojo is a degenerate band edge, then Vg{oJo) —>■ 
oo as A ^ oo; 
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c) if V & C°° on (0, a) and on {a,L), 0 < a < L, and v{x) > c > 0 on one of these 
intervals, v(x) < —c < 0 on another one, then uq is not a degenerate band edge when A is 
large enough; 

d) if q has the form then oj = ojn = X’ non-degenerate band edges, and 

Vg{uJo) = 0{A~^) when A —>• cxo and A\ sina;Z| < 1. 


Proof. Let the assumptions of part a) hold. Then there are two solutions = 'ipg, 
s = 1, 2, of the equation 

—'ll)” + {Av{x) — uP‘)i> = 0 (35) 

which have the following asymptotic behavior as a; G [0, L], |a; — a;o| < 1, A —>• cxo (see for 


example, P or [TT)]L 



fjiix) 

00 

~ e^o^AA)dtJ2aj{x,uj)A-^/^, ao = v-^/\x), 

0 


ij2{x) 

CXD 

~ e“fo a;)H“-^/^, /3 q = u“^/^(a:). 

(36) 


0 


This allows us to write the asymptotic behavior of the Priiffer monodromy matrix and 
its trace 2F{uj). We get that 

+ [^1'-") V^)dt){l + 0(A-'/^)). (37) 

Hence, |T"(a;)| > 1 when \uj — cuol < 1, and A is large enough. Thus, ujq belongs to a gap. 

Let the assumptions of part b) hold. Formulas (pF)|l still provide asymptotics of two 
solutions of (ESI), which oscillate as H > cxo, since now a/u is a purely imaginary function. 
This leads to the following analog of EH): 


F{ix) = (I 


t^(O) 

v{L) 


1/4 


+ 


AL) 

'u(O) 


|^/'^)cos(/ aJA\ v{f)\dt) + 0{A ^/^)), 


(38) 


and this formula can be differentiated with respect to uj. The main term of the asymptotic 
expansion of F{u}) does not depend on uj, and F''{uj) —>• 0 as H > cx). Together with Theorem 
Elthis proves that Vg{ujQ) —>• cxd as H —>• cx). 

Let the assumptions of part c) hold, and let u > c > 0 on (0, a), u < —c < 0 on (a, L) 
(the other case can be considered absolutely similarly. We shall use solutions 'ifiA 2 order 
to construct the transfer matrix over the interval (0,a), and similar solutions, with the 
lower limits in the integrals replaced by a, in order to construct the transfer matri^T^ over 
the interval {a,L). Then we arrive to the following formula for the main term of the 
asymptotic expansion of as H —»• cxo : = T^T^, where 


A 

rjn2 


^ [u(0)u(a_)]i/^ sinh Q cosh Q 

1^/^ sinP 


t (°+) 11/4 p 

I v(L) I cosr^ 

— — \v{a+)v{L) I sin P 


X 4 I v{a+)v{L) 
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Here 


v{a±) = lim v{x), Q 

x^azkO 



\/Av{t)dt, P 



\/ A\v{t)\dt. 


If a; = a;o is a degenerate band edge, then and ha^ to be equal to the identity 
matrix (Lemma Q). We evaluate non-diagonal elements of using formulas above, and 
equate them to zero. Then we arrive to a homogeneous system of equations for sinh Q cos P 
and cosh Q sin P with a non-zero determinant. Hence, 


sinh Q cos P = cosh Q sin P = 0. (39) 

Since Q > 0, we have sinhQ 7 ^ 0, coshQ 7 ^ 0, and therefore, can not be valid, i.e. our 
assumption on the existence of a degenerate band edge is wrong. 

Now we are going to prove the last statement of the Theorem. Let the potential be given 
by (jHl. Then the Priiffer monodromy matrix is the product: 

= T^(/+, L-)T^(r, P)T^(0+, r)T^(0-, 0+), 


where T^{a,b) is the transfer matrix which maps the Priiffer Cauchy data at a: = a of any 
solution to the equation i^ifo the Priiffer Cauchy data of the same solution at 

X = b, and the upper index ± in the arguments of indicates the limit value of from 
the right or the left, respectively. Elementary calculations give 


cos ul sin ujl 


■ 1 O' 


cos a;/ sin a;/ 


1 

0 

_ 1 

— sin ujl cos Lol 


^ 1 

L u) J 


— sina;Z cos a;/ 


1 

L oj J 



(40) 


and 

1 

F{oj) = -TrM^ = cosojL — —^sin^o;/. (41) 

One can easily check that is the identity matrix at a; = a;^ = n > 1 . Hence, all these 
points ujn are degenerate band edges. If H —>• cxo and A\ sina;Z| < 1 then (IHn|l and m imply 


|L sinwL^ sino;/coscnZ — ^ sin^ a;Z| A\A -|-coscuZ — sina;Z| | 

\k (o;)| =- 1 ^ ^ = - / — = 0[A). 

y sin^ ujL + ^ sin^ ul — sin"^ ujl ^ y o;^ -f 1 — ^ sin^ ujl 


This justihes the last statement of Theorem and completes the proof of the Theorem. 

Theorem El together with the last statement of Theorem El can produce an impression 
that both the group velocity and the reflection coefficient are small for the potential dSl 
if the spectrum of the incident wave packet belongs to a small enough neighborhood of a 
degenerate band edge a; = a;„ = ^,n>l. This impression is wrong. Let us consider this 
situation in more detail, taking into account that the wave packet does not have a hxed 
frequency a; = a;„, but the frequencies of the waves in the wave packet belong to an interval 
A centered at a;„. We assume that A belongs to the union of two bands adjacent at a; = a;„. 
Since |P(a;)| < 1 on these bands, (14111 implies that the length of these two bands has order 
0{A~^). We assume that A covers only a small part of the bands, since the wave packets 
with frequency spectrum strictly inside of a band will be studied below. So, we assume that 
the length |A| of A has order o{A~^). 
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If the spectrum of a wave packet belongs to an interval of length |A|, then the majority 
of the energy of the packet at a hxed time {t = 0) is concentrated on a space interval of 
length 1/|A|. The incident wave packet comes from a; = — cx), it propagates with the speed 
one in the region where the potential is zero (along the negative semiaxis, in particular), and 
it reaches the point x = 0 at time f = 0. After that, time 1/|A| is needed for the incident 
packet to enter the interval [0, NL] where the potential is supported. We would like to 
the time of the delay of propagation of the packet through the interval [0, NL] to be at 
of the same order. This leads to the relation i our case it means that 

'"“"O' 

Thus, we shall assume that |A| = o{A~^) and N satishes 

Theorem 7 Let the potential q have the form where A —>• cxo, let u E A where 

an interval centerd at ui = uin = O’ 01 = Then 

4 

"" uj-^A^sin^Nk (a;) +4’ 

where cos k{ijj) = F (a;) with F (a;) given in (O- 

Remark. One can see that the complete transparency (|tAr (a;) | = 1) takes place when 
sin Nk (cn) = 0. There is a hnite number (of order o{N)) of such points u in A, and to = to^ 
is one of them. Some very small neighborhoods of these points are also transparent enough. 
However, |tjv (i^) I ~ 1/A on the major part of A. This indicates that the majority of the 
energy of the incident pulse will be reflected. Perhaps, a more rigorous analysis of the 
situation is needed to justify the last statement. We shall provide this analysis below in a 
similar situation of the wave packets with the frequency spectrum in a narrow band. In this 
subsection, we shall restrict ourselves to considering the case of N = oo. 

Proof. One can easily evaluate the product (HOI) and then hnd ||M^|p - 2. We get 

- 2 = (— H--)sm^a;/. 

Since a; G A, we have that A sin a;/ —0. Using (ED we get 

A2 

sin^ k (u) = 1 — (a;) ~ — sin^ ul as A —oo, A sina;Z —> 0. (43) 

The statement of the Theorem follows from (EID and the last two formulas. 

We shall conclude this subsection by considering the limiting case of A = oo (the potential 
is supported on semiaxis x > 0), which supports the conclusion above that the majority of 
the energy of the incident wave packet is reflected. 

Theorem 8 Let 

j2 °° 

^ W 2 + nl) — 5{x — 2n/)], A —> oo. 

Lixib 


have 

least 


(42) 


A is 
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Let A be an interval centered at ujn = ^ and A| —0. Then the following relation is valid 
for the reflection coefficient: 

r{uj) = —1 + A —>• cx), a; G A. 


Proof. Formula (j24jl implies 


r(a;) + 1 


2 sin k (o;) + 2P — i {a — 6) 
2 sin k (a;) + (/? — 7 ) 


One can get easily from iOD that 


(44) 


A A A"^ 

a — 5 = —2—sina;Z(l + 0 ( 1 )), fl = o{—sinul), /? — 7 = — sina;/(l + o(l)). (45) 

OJ UJ 

The statement of the Theorem follows immediately from dH, (jlSl), and dH. 

Wave packets with frequencies inside of a band. For the same reasons as in the case of 
degenerate band edges, it is natural to consider the class of potentials of the form q{x) = 
Av{x), where n is a piece-vise smooth function, A —»• cx). If v{x) > cq > 0 then (see Theorem 
ini) any hxed cu = cjq belongs to a gap when A is large enough. This case is studied above (see 
Theorem n. If v{x) < —Co < 0 and v is smooth enough then F{uj) is given by Thus, 
|F(a;)| can be made less than one when A is large, and the chosen a; will belong to a band. 
However, dHHll implies that Fflco) = 0{A~^^‘^) as A —>• cxo. Since coskfjj) = F{uj), we have 

k'{Lj) sin k{Lj) = 0{A~^^'^), A —>• cxo. (46) 


If uj is strictly inside of a band then sinfc(a;) is separated from zero, and (I46j) implies that 
k'{u) —>• 0 as A —»• cx), and therefore, the group velocity tends to inhnity. 

A natural way to achieve a small group velocity is to consider a medium with a narrow 
band. Since k{uj) changes by tt over a band, k'{uj) has an order 0{7i/d) in a band of width 
d. Hence, the group velocity will have order 0{Ld/7r) and it is small when d is small. The 
easiest way to get a narrow band is to consider a potential q{x) = Av{x), where n = 0 on 
some interval (a, b) C (0, L) and v{x) < — cq < 0 on the remaining part of the interval of 
periodicity. Solutions ^ G of the equation 

d^ 

Htlj — Xf) = f, where H = — — -|- Avix), 1mA 7 ^ 0, / G 

dx^ 

tend to zero on the set {x : v{x) 7 ^ 0} as A —>• cx). This leads to the fact that the spectrum of 
the operator H when A is large enough, is close to the spectrum of the operator LLg = 
on the space Lf{U) with the Dirichlet boundary conditions on the boundary of U, where the 
set U consists of the interval (a, b) and all its shifts by a multiple of L. Thus, the spectrum 
of Hq is the set of eigenvalues of the Sturm-Liouville problem on (a, b) of inhnite multiplicity. 
Since the spectrum of H has a band-gap structure and the spectrum is close to the spectrum 
of FTo; the bands of H are narrow when A is large enough. 

Hence, we have a class of potentials for which the operator H has narrow bands, and 
therefore, the group velocity for frequencies inside these narrow bands is small. It remains 
to find out if long but finite slabs of such media are transparent enough (i. e. the reflection 
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is not too large). Unfortunately, the answer to the question above is negative: the smaller 
group velocity is in these narrow bands, the larger reflection is. This situation is similar to 
what we saw for the potential (El near a degenerate band edge. 

In order to make arguments simpler we restrict ourselves to considering a model case 
when the Hamiltonian has the form 

+ (47) 


which was studied in E). Note, that the interval (a, 6) in this case coincides with (0,L). 
Similarly to dUl) we get 




cosctiL sincuL 


1 

o 


cos ujL + - sin ujL sin ujL 

Ll) 

— sin ojL cos ujL 


^ 1 
La; J 


— cos ujL — sin ujL cos ujL 

L u) J 



and 


1 A 

F{(jj) = -TrM^^ = cosojL H-sin a;/. 

2 261^ 


(49) 


Theorem 9 1). If L and n are fixed, A —»• oo, then bands bn for the Hamiltonian have 
the form: 

T? TT 4 T? TT 

bn = [ 0 J 2 n-l, 0 J 2 n] = [-^{1 — + 0{A ^)), —], n=l,2... . 

2) . The group velocity Vg{uj) = L/k'{u}) on each band bn has the estimate 

IWI< ^(1 + 0(44-^)), A^oo. (50) 

3) . The absolute value of the transmission coefficient for iv E bn is equal to 


|tAr(a;)| 


1 



sin Nk{Lj) 1 2 
sin k(Lo) ' 


5 


and it has the order 0(A ^) for all iv E bn, n > 1, except for very narrow (of order 0(A ^)) 
neighborhoods of N — 1 transparency points where 


smNk{u)=0, sin/c(a;) 7 ^ 0. 


Proof. The inequality |T’(a;)| < 1 can be valid only if ujL is close to nvr, n > 1 (since 
a;“^sina;L has to be small). In order to hnd the bands, we put uiL = nvr + a, |cr| << 1. 
Then cosujL = (—1)"cos(T, sina;L = (—1)" sin cr, and the inequality |F| < 1 takes the form 


|^(c^)|<l, kl « 1, 


where g{a) 


COSCT + 


AL 

2{n'K + cr) 


sin cr. 


(51) 
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Note that ^'(O) = 1 and g'icr) > cqA > 0 when |(t| < vr/d and A oo. From here it 
follows that m holds if —e > cr > 0 where £ is the solntion of the eqnation g{e) = —1, i.e. 

AL 2 / 

0 = l + cos£- ;--Sin£ = 2 COS e 2 

2{n7i - e) ' 

, , , AL sin e/2, 

= cose/2 2cose/2-. 

rni — e 

Thns e = e„ is the solntion of the eqnation 

tane/2 = (52) 

i.e. En = ^ + 0{A~‘^). Since bn = [(nvr — en)/L, mi/L], the hrst statement of Theorem IHl is 
proved. 

Let ns prove the second statement. One can easily check that 

A T? 

F'Ujj) = {—lY -(1 + O(A“0) when a; G 6n, A —>• cx). (53) 

2mi 

Then from m it follows that 

A T? 

/c'(a;) sin/c(a;) = (—-(1 + 0(A“^)) when a; G A —>• cx), (54) 

2mr 

which implies the second statement of the Theorem. The last statement follows immediately 
from (ESI) and (gOl). 

Theorem n is proved. 

Theorem El implies that if the freqnency spectrnm of an incident pnlse belongs to a band 
bn of the operator (gZl) then only a small part (of order 0(A ^)) of the energy of the pnlse will 
propagate throngh the medinm, and the main part will be reflected. This statement needs 
to be justihed more rigorously, in particular, because of the existence of the transparency 
zones. Besides, the question about the speed of propagation of the transmitted part of the 
energy remains to be addressed, since the usual arguments about the group velocity concern 
the propagation of wave packets in inhnite media, not through a hnite slab of a periodic 
medium. The next theorem deals with both problems. In particular, we shall show that the 
main part of the incident pulse will be reflected. Only a part of the energy of the pulse of 
order 0(A“^) enters the media. It propagates with the group velocity, and only a negligible 
part of order 0(A““) moves faster. 

The following non-stationary problem describes the propagation of the pulse u = uo{x—t) 
through a hnite slab of the periodic medium: 


AL sine/2 cose/2 
nn — e 


N-l 

Utt = Uxx — A''^^ 6 {x — nL) u, {t,x) E u = Uo{x — t), t<0, (55) 

n=0 

where 

uo{x) = B-^/^a{x/B)R^^°, a G C'o“([-l, 0]), 5 = 0(A^+"), e > 0, (56) 
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and ujq belongs to bn- To be more exact, we assnme that 


TITT 0 

a;o = ^(l-^), O<0<4, n>l, (57) 

where 6 does not depend on A. The form of Uq and the meaning of the restrictions dsni), 
(EZI will be explained below, bnt hrst let ns mention that the energy conservation law holds 
for the solntions u of the energy 





N-l 

+ \ut\‘^)dx + A \u(mL, 

m=0 


does not depend on time. Since the energy of the solntion at t 


0 is eqnal to 


E., 


no 



—Uo{x)\‘^dx, 

dx 


(58) 


we have 

E = E 

for all f > 0. 

The fnnction uq is normalized in snch a way that the energy of the incident wave UQ{x — t) 
at f = 0 is separated from zero and inhnity and tends to a constant as B —>■ oo : 


Euo =‘2ujl\\a{x)\\l2 + 0{B ^). (59) 

The Fonrier transform of Uq is eqnal to 

Uo{uj) = B^^‘^a{B{uj — cuo)). ( 60 ) 


From Theorem M and (EZD it follows that Uq belongs to bn and is located not very close to 
the end points of bn- The fnnction uq is snpported, mostly, on a small interval 

A = {a; : \uj — cuqI < (61) 

located inside of bn- In fact, from Theorem El and (I57|l it follows that 

A C bn and p{A.,dbn) > fd/2A ii A > { 2 L/ ^ (^52) 


where fd = min(0, A — 9) and p(A, dbn) is the distance between A and the end points of 
bn- Fnrthermore, since a G (7“, we have 5(a;) = 0(|a;|“°°) as |a;| —>• oo- This and the 
assnmption B = 0{A^^‘^) imply the following estimate, which holds when a; ^ A : 

A l+e/2 

|S„(a,)| = |B‘-'=5(B(a,-i.j„))| = B‘-'"0((—A ^ oo. (63) 

Hence, Uq can be considered as a wave packet with freqnencies from bn- Let ns mention that 
it is impossible to choose snch a Uq that Uq is completely snpported on bn, since Uq would be 
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analytic in x in this case and conld not be eqnal to zero for x > 0, while the incident wave 
Uo{x — t) is supported on the semi-axis x < 0 when t < 0. 

We shall assume that N is not very small, since otherwise it would be unnatural to discuss 
a delay in the propagation of the incident pulse of space length B = > oo. One 

could consider values of N which are independent of A and B, but we shall assume that 
N = 0{A'^) with some 'j > e. Then the hnal estimates will not depend on N, so the result 
will be simpler. 

Let us evaluate the group velocity Vg{uj) = L/k'{(jj) for wave packets with frequencies 
from A. From (gni), dszD and p[|) it follows that 

cos k{uj) = - ^) + when a; e A. 

Hence, 

I sin k{!jo)\ = -\/A9 — 6^ + 0{A~'^^‘^) when o; G A. (64) 

From here and (1^ it follows that 




AL 


Thus, the time 


t = to(l + 0{A where to = 


ANB^ 


n'K\/ A6 — 

is needed for a signal to pass through the interval [0, AL] of the periodic medium, if the 
signal propagates with the group velocity V^(a;o). 


Theorem 10 (see m) The following estimates are valid for the solution u to the problem 
(Eil when t < fo(l — d), 6 > 0, and A —> oo: 

1) u = uq{x — t) — uo{—x — t) + g{x + t) when x < 0, where 

/ OO 

\g\x)\^dx < 6 * 2 ^“^ < cxo. 

■OO 


OO 

[ (ul + u‘()dx = 0{A-^) 

NL 
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